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A3S TRACT 

This paper preseats a statistical software package 
developed for use on the Apple II Plus microcomputer, 
modified with the AppLe Pascal language card. The program 
addresses the following: determination of confidence 

intervals for single and bivariate populations; hypothesis 
testing for one and two parameters; computation of 
cumulative distribution values for the Normal, Student's T, 
Chi-square, F, Binomial and Poisson distributions; 
computation of quantile values for the Normal, Student's T, 
Chi-square and F distributions. The program aLso has the 
capability to store, retrieve, and modify data for use with 
the statistical procedures. The program was written in (JCSD 
Pascal, which because of its portability indicates that 
little or no modification would be required to use it with 
other computers which are UCSD Pascal compatible. In 
addition, because of Pascal's block structure, the program 
can be easily modified or enhanced to include other 
statistical procedures which are of interest to the user. 
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I. 



INTRODUCTION 



For the operations analyst or for that matter anyone who 
utilizes the methodology and problem solving approach of the 
operations analyst, some form of computational device is a 
necessity. Cost, rigid interface requirements, and a number 
of other factors have in the past frustrated and stifled the 
analyst in bringing to bear the necessary computational 
power on his problem. In the late 1970's, however, Texas 
Instruments introduced a significant amount of computing 
power packaged in the programmable TI-59 calculator. 
Revolutionary is perhaps a bit too strong as a description 
of the impact that this and like devices have made in the 
operations research community, but certainly most analysts 
would agree that the amount of computing power that can now 
be held in one' s hand and used to solve problems is indeed 
remarkable. 

In spite of this impact, however, the hand-held 
calculator's contribution might even now be waning and 
yielding to a more spectacular capability found in the 
microcomputer. The growth in the capability of these 
devices since the beginning of the industry in 1971 has been 
phenomenal. tiastrakas details this growth and gives a 
glimpse of the possible direction that the 1980 's will see 
in this industry [Ref. ij . One factor accounting for the 
growth in microcomputer technology and one that serves to 
insure its future is the intense competition that pervades 
the industry. When one thinks of hand-held calculators, he 
thinks of Texas Instruments or perhaps Hewlett-Packard. When 
one thinks of microcomputers, he may think of \pple, Exidy, 
North Star, PET, T3S-30, or a number of other devices with 
like capability. 
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A typical microcoraput er package costing between 32,000 
and $2,500 might consist of the following: the computer 

(about the size of a small portable typewriter) with 64 K 
bytes of usable random access memory, two floppy-disk drives 
(5 1/4 inch diameter) for program storage, and a black and 
white monitor for output display. A variety of programming 
languages are also available including the University of 
California, San Diego ((JCSD) Pascal, BASIC, FORTRAN, PILOT, 
etc. Some of the microcomputers also have a graphics 
capability allowing the user to visualize mathematical 
forms, plot graphs, and plot data observations. 

The American populace is being conditioned through 
current periodicals and news features to expect an 
increasingly important role for the microcomputer in 
everyday life — from grocery shopping to environmental 
control for the home. The December 1, 1980 issue of "u. S. 
News and World Report” predicts that eighty percent of the 
United States' households will have a microcomputer by 1990. 
In view of the prospective proliferation of these devices 
and also the significant computing power that now exists, 
the operations research analyst can ill afford not to begin 
to exploit the capabilities of microcomputers. Indeed, the 
hardware capabilities have grown and are growing so rapidly 
that today good compatible software to support these 
impressive capabilities is seriously lacking in specialized 
fields such as operations research. It is this software 
deficiency that this thesis addresses. 
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II. BACKGROOND 



The operations analyst can, on occasion, be required to 
establish confidence intervals or test hypotheses about an 
unknown parameter from a known or assumed population. The 
computer program written as a part of this thesis allows the 
analyst to quickly and easily accomplish these tasks when 
observations are from Normal, Exponential or Bernoulli 
populations. 

A. DESCRIPTION OF THE MICROCOMPUTER SYSTEM 

The software development was done on an Apple II Plus 
microcomputer with two floppy-disk drives (5 1/4 inch 
diameter) as add-on peripherals. The system is equipped 
with a language system giving the capability to use the 
University of California, San Diego (UCSD) Pascal 
programming language as well as the BASIC language which 
comes resident with the computer. The standard Apple 
computer has a forty column output display which makes it 
compatible for use with a standard television set. This 
capability can be enhanced to an eighty column display with 
an additional peripheral device, provided a monitor is used 
for display in lieu of a television. The output format for 
the program is written for an eighty column display device; 
however, using the special built-in features of the Apple, 
the format can easily be made to display split screen on a 
system not equipped with an expanded display peripheral. 

B. SOFTWARE DEVELOPMENT 

Ling and Muller give several considerations which should 
be observed in the development of software for statistical 
analysis ^Ref. 2, Ref. 3j . Among these considerations are 
the following: 
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1 . Choice of Programming Language 

The Pascal programming language used in this 
software development effort has a number of features which 
make it particularly attractive for use with microcomputers. 
First, Pascal is a very concise language. The compiler is 
small and compact and fits easily within the available off- 
line storage space of the floppy-disk. Second, Pascal is a 
high-level, general-purpose language ^Ref. 'uj . The 
language was originally introduced in 1971, which is recent 
in comparison to most high-level languages. It was intended 
to be used in teaching new programmers good techniques and 
style. The Pascal language fully exploits the fundamental 
concepts of structured programming, which is a technique 
used by many professional programmers to write large complex 
computer programs ^Ref. 4 . The use of these techniques 
facilitates developing programs in a modular fashion (i. e. , 
break the overall package into logical sub-packages and 
proceed to program, debug, and validate each sub-package 
individually) . The final step in the process is to combine 
the sub-packages to form the overall package. Using the 
attendant statistical package as an example, six sub- 
packages or modules comprise the overall package. Each 
module is independent of the others and can stand alone when 
compiled with the main program. Third, the Pascal language 
performs arithmetic computations significantly faster than 
the BASIC language. Pascal is often implemented as a 
"pseudo" interpreted Language meaning that the text versions 
of programs are first compiled into a code file. It is 
during the compilation phase that syntax errors are detected 
by the compiler and brought to the programmer's attention. 
This code file is interpreted and executed during the 
execution phase of the program. Host processors can and 
typically do interpret a Pascal code file significantly 
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faster than a corresponding 3ASI3 program which performs the 
same computations ^Ref . 4j . 

2. Computational Efficiency 

In spite of thair impressive capabilities, 
microcomputers are decidedly inferior to largar computers in 
the two key areas of computational speed and accuracy. 
Pascal, as implemented on the Apple, will only perform 
computations using six decimal places of accuracy ^Ref. 5 ] 
and displays only five places past the decimal. Accordingly, 
algorithms which are used and work very well on larger 
computers might have no chance of producing the same results 
on microcomputers simply because the execution time is 
excessive, or they require double precision arithmetic. 

C. PORTABILITY 

A portable program is one that can be run on a number of 
different computer systems [^Ref. 2j . Programs written in 
the Pascal language are valuable in this regard in that they 
may be run on a variety of microcomputers without 
alteration. Some of the machines for which this is true are 
computers with the following microprocessors: 8080, 3085, 

Z80, 6502 (Apple), 6800, and 9900. Portability is a very 
strong asset of the Pascal language even though it is 
accomplished at the expense of reduced computational speed. 
The compiled cods version of programs is called ’'p-code”. 
Each machine has a special interpreter program which takes 

ti n , 

the p-code and converts it to a form compatible with the 
existing host microprocessor. 

D. EASILY RETRIEVABLE INFORMATION FILES 

Ling suggests the use of help files to allow the user to 
make efficient use of a program ^Ref. 2j . The number of 
information files which are available and their content is 
clearly a matter of judgment based on assumptions concerning 
the knowledge of the user population about the procedures 



used in the program. Anticipating that most of the users of 
the accompanying statistical package will be familiar with 
the basic concepts of the procedures themselves, only one 
help file which pertains to data entry requirements is 
included. 

E. SIMPLE USER INTERFACE 

The user is called upon to make numerical entries 
throughout the program. Because Pascal is a strongly typed 
language, a variable of type integer cannot be assigned 
floating point values. Indeed, this requirement is so 
strict that if the program expects the user to enter an 
integer and he accidently enters a number in decimal point 
notation, the program will abort and cause the entire system 
to re-initializ e. Situations such as this are, of course, 
undesirable. Moreover, all programs written with the intent 
of establishing a dialogue between the program and the user 
should be as trouble-free as possible for the user and 
minimize as much as possible the user's chances of 
committing a fatal error when responding to program prompts 
or entering data. In general, schemes to accomplish this 
are costly in terms of computational efficiency and 
programming steps; however, user convenience is almost 
always worth the costs. The problem with data entry or 
numerical entry from the keyboard is addressed in this 
statistical package by making all numerical entries using 
a string and converting the string to the numerical value 
it represents. A string is a variable type which is a 
linear array of characters. For example, given a string 
variable called 'S', which is assigned a value of *235', 
the characte 
•3', and 's[ 

'9' have corresponding numerical American Standard Code 
for Information Interchange (ASCII) values of forty-eight 



r in array position ' S^lj' is '2', 'S^j' is 
3 J 1 is '5'. Since the characters 'O' through 



through fifty-seven, conversion is accomplished by 
subtracting forty-eight from the character's ASCII value. 

A Pascal procedure which converts strings to numbers is 
shown in Appendix B. 

The advantages in using a scheme such as this for all 
numerical entries are two-fold. First, the user can take 
advantage of the direct cursor addressing available on many 
microcomputer video displays to correct a data entry prior 
to entering it into computer memory. This is not possible 
when the program expects real numbers as input. Second, the 
program segment which converts the string may return either 
an integer or a real number, whichever is required by the 
program. This avoids having the user concerned with the 
typing requirements demanded by the Pascal language. 
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HI. ALG0RIT8 MS 



Algorithms used for statistical computations on 
microcomputers should be selected with the goal of providing 
the best accuracy achievable at the minimum computation 
time. Given the limitation of only six decimal places of 
accuracy on the Apple microcomputer/ many algorithms 
requiring more precision in their computations must be 
rejected. The algorithms remaining as candidates must be 
carefully screened to insure that their required 
computations are reasonable from a time standpoint and that 
they exercise fully the accuracy capability of the machine. 



A. CALCULATION OF VARIANCE 



To illustrate 



accuracy 
[Ref. 6 



consider 



the complications stemming from reduced 
the following example given by Forsythe 



Find the variance of the following set of numbers: 
48499 , 48503 , 48500 , 48498 , 48500 
The definition for the variance is 




where S* = variance 

N = number of observations 
X = sample mean. 



This formula can be expanded to the following form: 



t 

S = 



1 



N - 1 





Use of this formula in the Apple, however, would yield a 
variance of zero, simply because the required accuracy is 
not available for computations. The answer of zero is, of 
course, incorrect. Focsythe in his article on statistical 
computing offers the following alternative algorithm to 
compute the variance. 



SUM : = 0.0; 

S2 := 0.0; 

READ (N) ; 

FOR I := 1 TO N DO 
BEGIN 

READ (X) ; 

DEVIATION ;= X - S(JM ; 

SUM := SUM ♦ DEVIATION/I; 

S2 := S2 + DEVIATION * (X - SUM) ; 

END; 

S2 ;= S2/(N - 1) ; 

a 

This algorithm produces the correct answer, variance = 
3.5 U . It also illustrates how many of the limitations of 
the microcomputer can be overcome through careful selection 
of algorithms. 

B. DISTRIBUTIONS AND INVERSES 

The algorithms used to compute probability distributions 
and inverses and the source of each are listed in Appendix 
A. When compared to the tabular values listed in Dixon and 
Massey |^Ref . 7 ] , the algorithms are accurate to at least 

three decimal places in probability with the exception of 
the F distribution. The F distribution is accurate to three 
decimal places in probability in almost all cases; however, 
some values may differ from the listed tabular values by as 
much as .002 in probability. Although the F quantiles may 
differ slightly from the listed tabular values, the 
probabilities corresponding to the values given by the 
algorithm are accurate to three decimal places. All of the 
algorithms produce results within three seconds except the T 
distribution. When computing for very large degrees of 
freedom for the T distribution, computation time is a 
function of the degrees of freedom. Typically 1000 degrees 
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of freedom takes approximately eignt seconds. 

Often, increased accuracy can only be gained through 
iterative techniques. This is the case with many of the 
probability distributions and their inverses, the 
T distribution being an example. In the algorithms shown in 
Appendix A much of the excessive computation time has been 
alleviated by combining two or more algorithms. For 
example, in the Chi-square distribution, the F distribution 
and the inverse F distribution, one algorithm is used for 
small degrees of freedom and another for large degrees of 
freedom. Good algorithms for these distributions exist for 
large degrees of freedom which are not based on iterative 
techniques and hence are computationally fast. However, for 
small degrees of freedom, their accuracy falls off rapidly. 
Conversely, the algorithms using iterative techniques are 
very accurate at all ranges, but slow for the larger degrees 
of freedom since the number of iterations required is 
proportional to degrees of freedom. The break between small 
and large degrees of freedom is purely subjective, based on 
choosing the best combination of speed and accuracy. 
Selection of the algorithms themselves was likewise based on 
the best combination of speed and accuracy. 

C. POISSON AND BINOMIAL DISTRIBUTIONS 

Large mainframes have the capability of computing the 
Poisson and Binomial Distributions directly from their 
definitions. 



POISSON: 

If the random variable X is distributed Poisson with 
parameter 

5 



then pf X <= kl = — — — 

L J 3 = o j! 
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BINOMIAL: 

If the random variable X is distributed Binomial with 
parameters [n,p j, 

p’ d-?r 

Both of these expressions contain factorials and 

summations; consequently for large values of 1 k’ in the 

Poisson and ’n' and "c 1 in the Binomial, execution time 

might be excessive, or the intermediate values in 

33 

computation might exceed ’10 1 (the maximum number capable 

of being represented oa the Apple) . A better solution for 
the calculation of these probabilites is to use the 
Chi-square identity for the Poisson and the F identity for 
the Binomial [aef. 8, Ref. 9j . 

POISSON: 

Given X ~ Poisson 

then Pr [x <= k] = 1 - X^k+z) (2 A) 

where 2k + 2 are the degrees of freedom 
of the Chi-sguare variate. 

BINOMIAL: 

Given X ^ Binomial [n, p] 

then P [x <= k] = P [P (2 „. u>il „ 2) <= ^ 

where 2n-2k equals the degrees of freedom for the numerator 
and 2k+2 equals the degrees of freedom for the denominator. 

When k = n, the probaoility is 1.0. 

D. STATISTICAL ALGORITHMS 

Derivation of the algorithms used in finding confidence 
intervals and hypothesis testing is shown in Storer 
[Ref. lo] . 



then pfx <= kl = 4 
L J j = o 
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17. DESCRIPTION OF PACKAGE 



The options available to the user are shown in the block 
diagram (Figure 1) . ihen the program is executed, the user 
begins in the “outer level”. To proceed he must select from 
one of the seven options which ace shown on a “menu" (Figure 
2 ) . 

A. THE DATA ENTRY MODOLE 

Selecting option '6' from the outer level menu (Figure 
2) will cause the data entry module menu (Figure 3) to 
appear on the screen. 

1 . General 

The data entry requirements are intended to be as 
trouble-free as possible for the user. The user is prompted 
for data input by the following line: 

Record N — > 

' N' is an integer sequentially updated by the 
program when the 'return' key is pressed. Following the 
arrow, the user inputs as many data values as he wishes 
with entries separated by one or more spaces. The only 
restriction is that he should not exceed the length of the 
display line. The nomenclature, record, indicates only a 
logical or convenient grouping of data from the user's 
point of view. 

Prior to entering any data, the user is asked 
whether or not the observations he is entering are paired. 
Since two of the statistical procedures are predicated on 
paired data, answering ' yes' to this question will cause 
summary statistics to be computed on 'X,Y' pairs. These 
summary statistics are only good for use with the 
procedures requiring paired observations. All data entered 
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Figure 1. Block Diagram 
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1) Hypothesis Testing One Parameter 

2) Hypothesis Testing Two Parameters 

3) Confidence Intervals Single Population 

4) Confidence Intervals Bivariate Populations 

5) Distributions and Inverses 

6) Data Entry 
Q) uit 

Figure 2. The Outer Level Menu 



1) Instructions 

2) Create a new data file 

3) Correct/add to existing datafile 

4) Enter data without storing 

5) Review existing datafile - 

6) Review summary statistics of existing file 
Q) uit and return to outer level 

Figure 3. The Data Entry Module Menu 
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must be in 'X,y pairs with the 'X' ooservation listed 
first. All other statistical procedures and their attendant 
data entry requirements assume that all of the observations 
come from a single population. These general requirements 
are available to the user when he selects option ' 1* from 
the data entry menu (Figure 3) . 

2. Create A New Data File 

Selecting option *2’ (Figure 3) will prompt the user 
to specify a file name for the data observations. After 
responding with a file name, the operating system 
establishes a directory entry at the beginning of the 
largest unused block of space containing at least fifteen 
blocks on the specified disk. Because the filing system in 
the University of California, San Diego (UCSD) Pascal 
language is random access, each data file entered will have 
allocated fifteen blocks of space to insure that there is 
enough room to extend the file if the user desires to do so 
at some later time. 

The floppy-disk (5 1/4 inch diameter) used by the 
Apple System provides a storage space of 280 blocks. This 
results in a capability to store seventeen data files on 
each disk. Bach data file can contain a maximum of ninety 
records. Since each record can contain as many observations 
as the user desires provided that it does not exceed the 
length of the display line, a reasonable planning figure is 
eight data observations per record. This results in an 
upper limit of approximately 720 data entries per data 
file. 

Data entry for a new file begins with 'Record 1'. 

The user's only concerns when entering data should be to 
separate each observation by one or more spaces, to not 
exceed the length of the display line, and to enter 'X,5f' 
pairs if he has previously indicated paired observations. 
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Since the data is entered as a string variable and then 
coaverted, it is possible to correct any entry prior to 
going to the next record by simply using the 'back arrow' 
key to move to the place at which the correction is to be 
male. The legal symbols which may be used in data entry are 
the digits 'O' through '9', the decimal point, the comma, 
the plus (+) and minus (-) signs, and the ' E' for scientific 
notation. Commas are ignored by the procedure which 
converts the strings and are included only as a convenience 
for the user. 

Pressing the 'return' key at the end of a record 
terminates that record and prompts the user to input the 
next record in sequence. If the user inadvertantly enters 
an illegal character while entering data, the program will 
advise him of this, indicate what the character was, and 
prompt him to reenter the record. 

Once the user has entered all of his data, he must 
press the 'escape' key and the 'return' key following the 
last data entry. This will close the data file. When the 
file is closed, the following summary statistics are 
computed on the observations in the file: 

a. The Sample Mean 

b. The Sample Standard Deviation (N-1) 

c. The Sum of the Observations 

d. The Sum of Squares of the Observations 

e. The Number of Observations 

For paired data, the same* statistics are computed; 
however, they are computed on the differences of the 'X, y 
pairs. Hence, the number of observations for a data fiLe of 
paired observations is exactly half of the total number of 
data entries. 
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The summary statistics for each data file are kept 
in a separate file that is initially established on the same 
disk as newly created data files. The file of summary 
statistics requires one block of disk storage space. It is 
differentiated from the original file of observations by an 
•S' concatenated to the original file name. 

3. Correct/&dd to Existing Data File 

Selecting option '3' (Figure 3) while in the data 
entry module will prompt the user to specify a file name. 
Once the data file is retrieved, the monitor will show the 
name of the file, the number of records in the file, and the 
number of observations. Immediately after this information 
is the following prompt: 

Enter Record (lumber — > 

Selecting any record number between • 1 • and the 
total number of records in the file will cause the retrieval 
of that record and will display as follows: 

old Record: 

Record 2 — > 12345 

You may replace the complete record or 
press <RTN> for no changes. 

Record 2 — > 

Pressing 'return' leaves the existing record 
unchanged and prompts the user to enter another record 
number. If corrections are to be made, the user must enter 
the complete new record opposite the lower prompt and press 
'return*. The updated record will display as follows: 

Record 2 — > 12335 

Press <RTN> if OK, <E5C> if not. 

Pressing 'return' completes the update and prompts 
for a new record. Pressing 'escape' produces the following 
display: 
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Enter corrected data and <RTN> 
Record N — > 



Ihe user then may retype the line. This sequence 
may be repeated until the record appears as the user wants. 
Entire records are erased or deleted by typing a space to 
produce a blank line. If the user desires only to update an 
existing data file and not extend it, he types the number, 
•-I*, in response to the prompt for record number, whica 
closes the updated file. 

If the user selects a record number that is greater 
than the number of records in the file, ha reenters the data 
entry phase beginning with the record number immediately 
following the last existing record. For axaraple, if there 
ace thirty-seven existing records in the file, selecting any 
number greater than thirty-seven will produce the following 
display: 

Enter new records. 

Enter <ESC> as the last entry and <RTN > to terminate input. 
Record 38 — > 

From this point, the user proceeds exactly as if he 
were in the data entry phase and terminates by pressing the 
•escape' key immediately following the last entry. 

4 . Enter Data Without Storing 

This option allows the user to compute the mean, 
standard deviation, sum of observations, sum of squares of 
observations, and number of observations. The format for 
entering data is exactly the same as previously discussed 
for creating a new file. As indicated, data is not stored 
under this option; hence, once a record is terminated, the 
entries for that record cannot be recovered nor changed. 

This option provides an expedient way to determine the 
summary statistics of a group of observations. 
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5 . Review an Existing Data Pile 

Selecting this option allows the user to quickly 
review any or all of an existing data file in blocks of ten 
records at a time. The user is prompted to enter the file 
name of an existing file. The program retrieves that file 
and displays the first ten records. Pressing 'return' at 
this point causes the next ten records to be displayed, 
etc. Pressing the 'escape' key at the end of any display, 
returns the user to the menu for the data entry module 
(Figure 3) . 

6. Review Summary Statistics of an Bxisting File 

All files that are stored on disk have an associated 
summary statistics file that is created by the program when 
a newly created or updated data file is closed. This file 
contains the following information: 

a. Hean 

b. Standard Deviation 

c. Sum of Observations 

d. Sum of Squares 

e. Humber of Observations 

The summary statistics file only is called when 
specifying a data file to be used in the other modules 
containing the statistical procedures. Because of this, it 
is not necessary to have the original data file on-line when 
performing the statistical procedures; only the summary 
statistics file is required. When using the filing system 
resident in the Apple Pascal language system to obtain 
directory listings of various disks, the summary statistics 
file is distinguished by a concatenated 'S' on the end of 
the original file name. For example, if the original file 
name was STATrDATAI , then the corresponding summary 
statistics file is named STAT:DATA1 S. Since the length of 
all summary statistics files is one block, it is possible to 
access a maximum of 27'4 files of summary statistics on any 
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one disk of storage. 

7. Q) nit 

Pressing the ' Q' key will return the user to the 
outer level (Figure 2). 

B. DISTRIBUTIONS AND QUANTILES MODULE 

Selection of any of the options from the distributions 
and guantiles module menu (Figure 4) will produce further 
prompts which require the user to enter the necessary 
information concerning values of the random variable, 
degrees of freedom, and probabilities as appropriate. 

After each computation the following prompt appears: 

C) ontinue or Q) uit 

Pressing the ' C* key will allow the programmer to 
continue calculation in the previously selected distribution 
or quantile. Pressing the ' Q' key will return the user to 
the menu for this module (Figure 4) . 

The Pascal language system allows the user to develop 
his own specialized libraries of often used subroutines for 
general-purpose or special-purpose computations. It is in 
such a library that the algorithms for the distributions and 
quantiles are kept. Using a special library has two major 
advantages. First, when the user is developing the main 
program, he is not penalized by extra compilation time for 
any of the routines in the library. The code in the library 
is linked by a separate process |^Ref . 4j . Second, the 
Libcary can be .used by other programs which require the use 
of the algocithms contained therein. The algorithms foe the 
distributions and quantiles fit logically in a library since 
it is likely that other statistical packages will require 
their use. The Apple reference manual explains the 
procedure used to establish new libraries Ref. 11 . 
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1 

2 

3 

4 

5 

6 



Normal Distribution 
T Distribution 
Chi-Square Distribution 
F Distribution 
Binomial Distribution 
Poisson Distribution 



A 

B 

C 

D 



Normal Quantiles 
T Quantiles 
Chi-Square Quantiles 
F Quantiles 



Q) uit and return to outer level 



Figure 4. The Distributions and Quantiles Module Menu 



POPULATION ASSUMPTIONS 



Population 



Assumptions 



1) Normal 

2) Normal 

3) Normal 
4} Normal 
5} Bernoulli 
6) Exponential 

Q) uit and return to outer level 



u?, sigma-sqr known 
u S sigma-sqr ? 
u & sigma-sqr ? 
u known, sigma-sqr ? 



u ==> mean 

sigma-sq ==> variance 



Parameter 



u 

u 

sigma-sq 

sigma-sq 




Figure 5. Confidence Intervals Single Population Menu 
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C. CONFIDENCE INTERVALS AND HYPOTHESIS TESTING MODULES 

Selecting either the confidence intervals or hypothesis 
testing options will cause the menu for that module to be 
displayed. All of these menus contain information similar 
to that shown in figure 5 for confidence intervals. 

The parameters about which an interval is to be computed 
or a hypothesis is to be made are listed in the right hand 
column of the menu. The distribution from which the 
observations came and assumptions about the populations are 
listed in columns one and two r respectively. 

1 . Data Requirements 

When one of the options is selected that does not 
involve Bernoulli or Poisson populations the following 
display appears: 

1) Use existing data file 

2) Enter data and store it 

3) Enter data w/o storing 

4) Use summary statistics 

Because of the nature of the observations, any tests 
or intervals involving Bernoulli or Poisson observations are 
entered using summary statistics; hence, for these cases, 
this display is skipped. 

a. Use existing Data File 

The user is prompted for the name of the data 
file. When the file name is entered, the program retrieves 
the summary statistics file associated with that file from 
disk. In the case of bivariate populations, two file names 
are needed, the first containing the 'X' observation and the 
second containing the • Y' observation. Paired observations, 
as noted previously, are entered in one file. 

b. Enter Data and Store it 

Selecting this option will display a short 
message to the user informing him that all data storage must 
be accomplished from the data entry module (Figure 3) . The 
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user has the option at this point to enter the data without 
stocing it or to return to the outer level (Figure 2) and 
select the data entry module. 

c. Enter Data W/0 Storing It 

Data is entered in the same format as 
previously discussed in the section pertaining to the data 
entry module. No permanent disk record is made of the 
entries. Hence, once each record of observations is 
terminated, there is no way to retrieve it to make 
corrections. Where bivariate populations are required, the 
user is prompted to enter all of the 'X' observations first 
and all of the 'X' observations second. For paired 
observations all of the 'X,X' pairs are entered as one 
population. 

d. Use Summary Statistics 

On all other tests or intervals, use of summary 
statistics is optional except as previously mentioned for 
data from Bernoulli or Poisson populations, at each 
prompt, the user is asked only for the necessary 
information to perform the statistical procedure he has 
selected. I he distinction is made in each prompt whether 
or not the statistics required are the sample parameter 
values (sample mean, sample standard deviation) or the true 
parameter values (true mean, true standard deviation) . 

2. Confidence Intervals 

When computing confidence intervals, the user must 
supply the additional information concerning the desired 
level of confidence and the type of interval (two-sided, 
one-sided upper, one-sided lower) . Typically, computation 
of the desired interval takas one or two seconds and is 
displayed as follows: 
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95.00 Percent confidence intervals for u 

Sample dean = 3.000 

Standard deviation = 1.581 

Upper = 4.964 

Lover = 1.036 

Another interval using same data, Y)es N) o — > 

Following the display, the user is asked whether or 
not he desires to compute another interval usiag the same 
data. If he responds ’yes', he may then vary the 
confidence level and/or the type of interval without 
having to again specify the data base. Answering ’no' will 
return the user to the menu for the module in which 
computations are currently being performed. 

3. Hypothesis Testing 

Hypothesis tests require the user to specify the 
null hypothesis. Typically, the hypotheses involve ’equal 
to', 'less than or equal to', or 'greater than or equal to' 
comparisons and are displayed for the user in a form 
similar to the one below: 

1) u = u [0l 

2) u <= u CO] 

3) u >= u CO] 

The symbol [o] represents the null hypothesis 
value. The user enters this value if required by the test. 
An example display following computation is as follows: 

HYPOTHESIS: U = U [0] 

Sample mean = 3. ODD 

u [Of = 2.500 

The P-value is 0.519 

Another test using the same data, Y)es N) o — > 

The user is not told to accept or reject the 
hypothesis; rather, he is given a p-value as shown above. 
The p-value, or probability level, is an indication of the 
level of confidence associated with the hypothesis ^Ref. 

12j. High p-values convey a high confidence in the null 
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hypothesis; conversely, low p-values reflect a lack of 
confidence in the validity of the hypothesis. A p-value of 
.05, for example, indicates that if the hypothesis is indeed 
true, there is only one chance in twenty that the data used 
in the test is consistent with the hypothesis. Upon 
completion of the test, the user has the option to perform 
another test with the same data or return to the menu for 
the module in which tests are currently being performed. 
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V. 



C0NCLUSI3NS 



The decade of the 1930's promises to be particularly 
bright in terms of affording the operations research analyst 
easy access to computing power. Microcomputers now 
available and their more capable descendants will doubtless 
play an important role in securing this access. However, 
equally as important are the software packages which will 
accompany these computers. 

The software package described in this paper provides 
the analyst with a useful set of statistical tools which can 
be used on one of the most popular, current microcomputers, 
the Apple. The five major modules contained in the package 
(confidence intervals for single and bivariate populations, 
hypothesis testing for one and two parameters, and the 
distributions and quantiles) are designed to be easy for the 
analyst to use and to cushion, as much as possible, 
potential user mistakes. The algorithms used throughout the 
program were chosen on the basis of being compatible with 
the microcomputer with respect to size and computing 
precision and providing the best combination of speed and 
accuracy. 

Pascal, the programming language used, offers not only 
the advantages of portability and increased computational 
speed, but also flexibility. Pascal is flexible in that 
large complex programs are programmed in modular segments 
which are then combined into the overall program. It 
follows that programs developed in this way are easily 
enhanced by the addition of new modules. Such is the case 
with this statistical package which could be significantly 
enhanced by the addition of a regression and aa analysis 
of variance module. 



This package and those to follow which ace compatible 
with current and future microcomputers can have a 
significant impact in the analyst community in two key 
areas. First, the analyst can be better educated. Simply 
alleviating the tedium which accompanies the application of 
many statistical procedures will give fledgling analysts the 
opportunity to work more problems and be exposed to a 
greater variety of situations in the school environment. 
Perhaps of equal importance, the educational environment can 
provide the opportunity to accustom the analyst to the 
capabilities that can and should be available for his use in 
a working environment. Second, by expanding computing power 
into areas which were not privileged to have it before, the 
educational process has a better opportunity to continue. 
Today, it is reasonable to assume that the professional 
growth of many analysts is stifled from a lack of computing 
machinery with which to attack his problems. 

Taking full advantage of the microcomputer's hardware 
capabilities requires efficient compatible software. In 
specialized areas such as operations research, the analysts 
themselves must logically provide the bulk of the effort in 
software development. The statistical package which is the 
subject of this thesis effort scarcely begins to provide the 
full complement of tools which the analyst requires. If the 
operations research community is to take advantage in a 
timely manner of the new opportunities afforded by 
microcomputers, effort must continue now in software 
development. 
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APPENDIX A 



FUNCTION Z ; 

/*##**##***##***# *# * # * * * * *** * ** ** **# * * 
** 

l'* 

i'* NORMAL DISTRIBUTION 

* SOME COMMON BASIC PROGRAMS 

■ '* 3RD ED., P. 128 

i’* 

l'* 

************************ ********** **** * 



CONST 
Cl = 
C2 s 
C3 = 
C4 : 
C5 = 



0.4361836: 

- 0 . 1 20 1 6 76 ; 
0.937298; 

0. 33267: 

2. 50662§3; 



VAR 

XX,XT,R,T : REAL; 



BEGIN 

IF STDEV <= 0.0 THEN 
BEGIN 
ERROR; 

EXIT (Z) ; 

END 

ELSE 

BEGIN 

XT * = X* 

X := ABS ( (X-MEAN) /STDEV) ; 
XX := x*x; 

R := EXP (-XX/2. 0) /C5 ; 

X := 1.0/(1 .0 + C4*ABS (X) ) 
T := 0.5 -S*(C1*X >C2*X*X 
IF XT < MEAN THEN 
Z := 0.5 - T 
ELSE 

Z : = T ♦ 0.5; 

END ; 

END; (* Z *) 



+C3*X*X*X) 
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###****»* 



FUNCTION INVZ 



************************ ******** ******** 



NORMAL QUANTILES 

HANDBOOK OF MATHEMATICAL FUNCTIONS 
P. 933 



**** ************** **** * * ********** *** * ** 



CONST 
Cl = 
C2 = 
C3 = 
D1 = 
D2 = 
D3 = 



2.515517; 
0.802853; 
0.010323 ; 
1. 432788; 
0. 189269 ; 
0. 001308; 



VAR 

PT,T,NUM,DEN 



REAL; 

OR (P <= 0.0) THEN 



3EGIN 

IF (P >= 1.0) 

BEGIN 
ERROR; 

EXIT (INVZ) ; 

END 

ELSE IF P > 0.5 THEN 
PT ;= 1.0 - P 
ELSE 

P T • — p • 

T := SORDINI 1 .0/(PT*PT) ) ) ; 

NUM := Cl + C2*T + C3*T*T; 

DEN ;= 1.0 ♦ D 1 *T + D2*T*T + D3*T*T*T; 
IF P > 0.5 THEN 

INVZ ;= T - NUM/DEN 
ELSE 

INVZ := - (T - NUM/DEN) ; 

END; (* INVZ *) 
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FUNCTION T 



**************************************** 



**** ************** ********************** 



CONST Cl = 0.63661977; 

VAR 

I,N : INTEGER: 

ANS,D2,F1,F2,T1 ,T2,XT : REAL; 
BEGIN 

IF DF < 1 THEN 
BEGIN 
ERROR; 

EXIT (T) ; 

END 

ELSE 

BEGIN 

XT := X; 



BEGIN 

ANS := 1. 0-Cl*ATAN (T1) ; 

IF DF <> 1 THEN 
BEGIN 

D2 := Z 1*T1 *T2 ; 

ANS := ANS - D2; 

END: 

IF DF <> 3 THEN 
BEGIN 

FI := 0.0: 

N := (DF-2) DC 7 2; 

FOR I := 1 TO N DO 
BEGIN 

F2 := 2 . 0*I-F 1 ; 

D2 := D2*T2*F2/ (F2+1) ; 
ANS := ANS-D2; 

END; 

END; 

END 

ELSE 

BEGIN 

D2 := T1*SQRT (T2) ; 

ANS := 1.0-D2; 

IF D2 <> 2 THEN 
BEGIN 

FI := 1.0: 

N := (DF-2) DIV 2: 

FOR I := 1 TO N DO 
BEGIN 

F2 := 2. 0*I-F 1 ; 

D2 := D2*T2*F2/ (F2+1) ; 
ANS := ANS-D2; 

END; (* END FOR *) 

END; 

END: (* END ELSE *) 

IF ANS <0.0 THEN 



* 

* 

* 

* 

* 

* 

* 



T DISTRIBUTION 
CACtl 

ALGORITHM 344 



* 

* 

* 

* 

* 

* 

* 
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AMS := 0.0; 

IP XT < 0.0 THEN 
T ;= ANS/2.0 
ELSE 

T := 1.0 - ANS/2.0; 

END; 

END; (* T *) 
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FUNCTION INVT ; 

/****************** ****** ****** ********* 



* 4c 

^ $ 

* T QUANTILES * 

* CACM * 

* ALGORITHM 396 * 

* * 

* * 



I**************** ** **** ** **** ****** **** * 

CONST 

HALFPI = 1.570796327; 

7AR 

PT,DEN,A,B,C,D,X, Y ; REAL; 



BEGIN 

PT * — p • 

IF (DF <: 1) OR (P >= 1.0) OR (P <= 0.0) THEN 
BEGIN 
ERROR; 

EXIT (INVT) ; 

END 

ELSE IF P > 0.5 THEN 
P := 2.0* (1 .0 - P) 

ELSE 

P ;= 2. 0*P; 

IF DF = 2 THEN 

INVT := SQRT(2.0/(P* (2.0 - P))-2.0) 

ELSE 
BEGIN 

IF DF = 1 THEN 
3EGIN 

P := P*HALFPI ; 

INVT := COS (?) /SIN (P) ; 

END 
ELSE 
BEGIN 
A 
B 
C 
D 
X 

y 

IF 
BEGIN 

X := IN vz (P*0. 5) ; 

Y * = X* X * 

I F* DF < $ THEN 

C := C + 0. 3* (DF - 4.5) * ( X + 0.6); 

C ;= ( ( (0. 05*D*X - 5 . 0) *i - 7 . 0 j * X - 2. 0^*X + 



= 1.0/ (DF - 0.5) ; 

= 48. 0/ (A* A) : 

= ((20700*A/B - 98.0) *A - 16.0)*A + 95.36; 

= ( (94 . 5/ (B + C) - 3.0) /3 + 1.0) *SQRT(A*HALFPI) *DF; 
= D*P; 

= XPN (X, (2 .0/DF) ) ; 

I > 0.05 + A THEN 



rJI< 

+ 1.0) *X; 



((0.4*? + 6. 3) *Y + 36.0) *Y + 



Y 

Y ;= A* (Y*Y) ; 

IF Y >0.002 THEN 

= EXP (Y) - 1.0 



B ♦ C; 
94.5) /C-Y-3. 0) / 



Y 

ELSE 

Y 

END 
ELSE 
Y : = 



IF PT >= 
INVT : = 



= 0. 5* (Y* Y) + Y; 



(1. 0/(((DF + 6.0) /(DF*Y)-0. 089*D-0. 8 2 2) * 
DF + 2. 0) *3.0) +0. 5/(DF + 4.0) ) *Y-1.0) * 
DF+1. 0^/(DF + 2.0) +1. 0/Y; 



5 THE 
SQRT (DF*Y) 
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ELSE 

INVT := - SQST (D?*Y) ; 
END; 

END; 

END; (* INVT *) 
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FUNCTION CHISQ ; 

r * ** *********** ** ** **** * * **************** 



* * 

$ # 

* CHI-SQUARE DISTRIBUTION * 

* SOME COMMON BASIC PROGRAMS * 

* (DF <= 40) P. 130 * 

* HANDBOOK OF MATHEMATICAL FUNCTIONS * 

* (DF > 40) P. 941 * 

* * 

* * 



******** ****** ******************** ****** 



VAR 

I : INTEGER: 

Y, POWER, TEMP, NU3, DEN, J,L,M 



< 1) OR (X <= 0.0) 



: REAL; 
THEN 



BEGIN 
IF (DF 
BEGIN 
ERROR; 

EXIT (CHISQ); 

END 

ELSE IF DF > 40 THEN 
BEGIN 

Y := ((XPN(X/DF, 1. 0/3.0) -1.0) + 2. 0/(9. 0*DF) | /SQRT(2.0) 

Y ;= Y*3. 0*SQRT (DF) ; 

IF Y < -4.3 THEN 

CHISQ := 1.0 
ELSE IF Y > 4.3 THEN 
CHISQ := 0.0 

ELSE CHISQ := Z (1,0. 0,1.0) ; 

END 

ELSE 

BEGIN 

DEN := 1.0; 

TEMP := DF; 

REPEAT 

DEN := DEN * TEMP: 

TEMP := TEMP - 2.0; 

UNTIL TEMP < 2. 0; 

POWER := (DF ♦ 1) DIV 2; 

NUM := XPN (X, POWER) *EXP (-X/2. 0) /DEN ; 

IF ODD (DF) THEN 

J := SQRT (2. 0/X/3. 1 415926) 

ELSE 

J : = 1.0; 

L := 1.0; 

M := 1.0; 

REPEAT 

DF := DF + 2; 

M := M* X/DF ; 

L := L+M: 

UNTIL M < 0. 0000001; 

L := L-M; 

CHISQ := J*NUH*L; 

END; 

END; (* CHISQ *) 
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FUNCTION IN7CHI 



**** ************** ********************** 
* * 

>jc 

* CHI-SQUARE QUANTILES * 

* CACM * 

* ALGORITHM 451 * 

* * 

* * 

* *************************************** 



CONSI 

Cl 


1 


1. 565326E-3: 


C2 


= 


1 .060438E-3; 


C3 


— 


-6.959356E-3 


C4 


— 


- 1.32329 3E-£; 


C5 




2 . 2775 79 E-2; 


C6 




-8. 986007E-3 


C7 




- 1 . 51390E-2; 


C8 


= 


2 . 5300 IE-3; 


C9 




-1 .450117E-3 


CIO 




5. 169654 E-3; 


C11 


= 


-1.15375 IE-2; 


C12 




1 . 1261 86 E-2; 


Cl 3 




2. 607083E-2; 
-8.42681 2E-4 ; 


Cl 4 




-2.23736 8E- 1 ; 


Cl 5 


= 


9 . 780499 E- 5 ; 
-8.553069E-3 


Cl 6 




Cl 7 




3. 125580E-3; 


Cl 8 


rr 


C19 




1 . 348028E-4; 


C20 


= 


4. 713941 E-1 : 
-1.425 296E-2 ; 


C21 


= 


1.0000886; 


A 1 




1 . 264461 E-2; 
-5 . 88609 OE-3 ; 


A2 


= 


A3 


= 


1 . 40C483E-2 ; 
-2. 304527E-2 


A4 




A 5 


25 


-1.091 121E-2; 


A6 




A7 


= 


3. 1354 1 1 E-3; 


A8 


5 : 


-2.728484E-4 ; 


A9 




-9.699681E-3 


A10 


= 


1. 316872E-2; 


A1 1 


= 


2. 6189 IE-2: 

3. 483789E-5; 
-8.729713E-3; 


A 1 2 




-2.2 22222E- 1 


A 1 3 




5. 40667E-5: 


A 1 4 




A15 




-7 . 27476 IE-4 


A 1 6 


= 


3. 2921 81 E-3; 


A17 




A18 


= 


4 . 7 1 4045E- 1 ; 


A 19 


= 


1.0; 















VAR 

F,F1,F2,TEMP : REAL; 

BEGIN 

IF (DF < 1) OR (P <= O.O) OR (P >= 1.0) THEN 
BEGIN 
ERROR; 

EXIT (INVCHI) ; 

END 

ELSE 

BEGIN 

P ;= 1.0 - P; 

IF DF = 1 THEN 
BEGIN 

TEMP ;= I NVZ (0.5*P) ; 

INVCHI ;= TEMP*TEMP; 

END 

ELSE IF DF = 2 THEN 
INVCHI := -2.0*LN(P) 

ELSE 

BEGIN 

F ;= DF; 

FI • = 1 0 /V • 

TEMP : = *INvi (1.0-P) ; 

F2 ;= SQRT(FI) *TEMP; 



IF DF 
TEMP 



ELSE 

TEMP 



INVCHI 

END; 



< ( 2 ♦ TRUNC (4.0*A8S (TEMP) ) ) THEN 

F2+C3) *F2+C4) *F2 

C7) *F1 + ( (((( (C8*C9*F2) *F2 
2+C12) *F2*C13) *F2 + C1 4) ) *FU 
♦F2+C 1 7) ♦ F2 + C 1 8) *F2 
♦ C21 






) *F2+C6) *F2 + < 
+ciO) *F2+c 1 1) *f: 
(Cfjf (C15 J *'F2 *C 161' 
+ C19) *F2+C2 0) *F2< 



:= ( ( (A1*A2*F2) *FH- ( 
♦ A5> *F2 + A6) ) ♦FI 
+ A10) *F2+A 1 1) *F 
A 1 4) ♦ F2+A15) *F 
+ A18) *F2+A1 9; 

:= TEMP*TEMP*TEMP*F; 



A3+A4 + F2) *F2 



if ((1 ( A7 + a 8*F2) *F2 +A9) *F2 
2* A12) ) *F1M( ( ( (A13*T 
2»- A16) *F2+A i 7) *F2*F2 



(* IMVCHI *) 



END; 

END; 
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FUNCTION FI; 

|' ** + * * ** ** *** ************ **** ************, 
'* ' * 

(4c * i 

(* F DISTRIBUTION * 



* SMALL DEGREES OF FREEDOM * 

* STATISTICAL COMPUTING * 

* P. 114 * 

* * 

★ *1 



(’**** *** ***** ** **** ****** ** * * ** ** * * **** * * 

CONST PI = 3. 14 159265; 

VAR A,A1 # AN1,AN2,ART,D1,D2,D3,R # S1,S2*DEL r XM,XK,T, 
TEM, C : REAL; 

I,M,N ; INTEGER; 



BEGIN 

AN1 := DF1 ; 

AN2 := DF2; 

A ;= DF1*X/(DF1*X+DF2) 
• = i - a* 

IF A 1 <= 0.6 THEN 
A 1 ;= 1.0E-37; 

D1 ;= AN1*0 . 5 ; 

D2 := AN2*0. 5 ; 

D3 ; = D 1 + D2 - 1 .0; 

R ; = 0.0: 

51 := 0.0; 

52 := 0.0: 

DEL := 1.0; 



0. 875 DO 



XM := 1.0; 

XK := 1.0; 

C := 0.25; 

N := DF2; 

HHILE C < 

BEGIN 

M := TRUNC ( D2) ; 

H := 2 * m; 

IF a = N THEN 
BEGIN 

N := TRDNC (D2) -1 ; 
IF N > 0 THEN 
BEGIN 
FOR I 
BEGIN 
SI 
D2 
D3 
TEM 



:= 1 TO N DO 
:= DEL+S1*R; 



= D2 - 1 
= D3 - 1. 
:= A1/D2; 
R := D3*TEM; 
S2 := (R+TEM> 
END; 

END; 

SI := DEL+S1*R; 
DEL := 0.0; 

T := -1.0: 

D3 := -1 . 6; 

:= A*S2; 

= C ♦ 0.5; 



0. 

0 ; 



*S2; 



S2 
C : 

END 
ELSE 
3EGIN 

N : = TRUNC jD2) ; 
IF N <> 0 THEN 



(* DEGREES OF FREEDOM ODD *) 
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3SGIN 
FOR I 
BEGIN 
SI 
D2 
D3 
TEM 



:= 1 TO N DO 



= DEL+S1 *H ; 

= D2 - 1.0; 

= D3 - 1.0; 

_ := A1/D2; 

R ;= D3 *T EM; 

S2 ;= (3 + TEH) *S2 
END; (* FOR *| 

END; 

51 ;= XK*S1 ; 

52 ;= XK*S2 ; 

ART ;= SORT If A 1) ; 

xa := xm* art; 

= (XM-ART) /A1 ; 

:= - 0 . 5 ; 

:= 2.0/PI; 

= C*2. 0; 



xa 

T 

D3 

XK 

C 

END: 
IF C 
BEGIN 
D2 
D3 
S2 
SI 
A1 
IF 



<= 0.875 THEN 



;= D1 ; 

:= D2 + D3; 

;= SI; 

;= 0 . 0 ; 

• r ^ • 

A 1 <i 0.0 THEN 
A 1 := 1.0E-37; 
;= DF1; 



N 

END; 

END; (* WHILE *) 

IF C < 1.125 THEN 

DEL := 4.0/PI*ATAN(T) : 
FI ;= XM*(S2 - SI) - DEL 
END; (* FI *) 



FUNCTION F2; 

(************ **************** ******************* 
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CONST 

Cl = 0 . 196854 ; C2 = 0.115194;C3 = 0.000344;C4 



VAR 

S,T,Z,Y,XX,YY,J,K 

BEGIN 



REAL; 



IF X > 
BEGIN 


1 THEN 


S : = 


DF1 ; 


T : = 


DF2 ; 


Z : = 
END 
ELSE 
BEGIN 


X; 


S : = 


DF2; 


T : = 


DF1 ; 


END; * 


1.0/X; 



J := 2.0/9.0/S; 

K ;= 2.0/9.0/T; 

Y := ABS( (1.0-K) *XPN(Z, (1. 0/3.0) > -1.0 + J)/ 
SQRT(K*XPN(Z, (2.0/3. 0) ) +J) ; 

YY := 1.0 + Y* (C1+Y* (C2+Y* (C3+Y*C4) ) ) ; 

XX := 0.5/XPN (YY,4.0) ; 

IF X >= 1.0 THEN 
F 2 j = 1.0 - XX 
ELSE 

F2 := XX; 

END; (* F2 *) 

FUNCTION F; 

BEGIN 

IF (DF1 < 1) OR (DF2 < 1) OR (X < 0.0) THEN 
BEGIN 
ERROR; 

EXIT (F) ; 

END; 

IF ( (DF1 < 100) AND ( DF2 < 100)) OR 
(DF1 < 20) DR (DF 2 < 20) THEN 
F := FI (X , DF 1 , DF2) 

ELSE 

F ;= F2 (X,DF1,DF2) ; 

END; (* F *) 
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CONST EPS = 0.005; 
TAR 



TEHP,PT, ENDR, ENDL,MIDPT,STEP, W,H,T,Y,Z 
1) OR (DF2 < 1) OR (P <= 0.0) 



REAL; 

(P >= 1 



OR (DF2 = 1) THEN 



BEGIN 

IF (DF 1 < 1) OR (DF2 < 1) OR (P <= 0.0) OR 
BEGIN 
ERROR; 

EXIT (INVF) ; 

END 

ELSE IF (DPI = 1) 

BEGIN 

IF DF 1 = 1 THEN 

INVF ;= SQR (INVT( (0.5*(1.0 + P) ) , DF2) ) 

ELSE 

INVF := 1 . 0/SQR (INVT ( ( 1 .O-P/2 . 0) , DF1) ) ; 

END 

ELSE 

BEGIN 

Z := INVZ(P): 

Y ;= (Z*Z - 3 . 0 ) /S . 0 ; 

H ;= 2.0 * 1. 0/ (1 .3/ (DF2-1. 0) + 1 . 0/ (DPI - 1 . 0) ) ; 

W := (Z*SQRT (H+Y) ) /H - ( 1 . 0/ (DP 1- 1 . 0) -1 . 0/ (DF2- 1 
jY + 5. 0/6.0 -2 . 0/ (3 . 0*H) ) ; 

T ;= EXP (2. 0*W) ; 

IF (DPI <= 15) OR (DF2 <= 15) THEN 
BEGIN 

STEP := 0 . 1 *T ; 

ENDR := T; 

REPEAT 

TEMP := ENDR; 

ENDL ;= ENDR - STEP; 

IF ENDL <= 0.0 THEN 
ENDL ;= 0.0; 

ENDR := ENDL; 

PT := FI (ENDL. DF1,DF2) ; 

UNTIL (PT <= P) OR ((TEMP - ENDL) <= EPS) ; 

ENDR := TEMP; 

IF ABS (ENDR - ENDL) > EPS THEN 
BEGIN 

WHILE (ENDR - ENDL) > 0.1 DO 
BEGIN 

T ;= (ENDR *■ ENDL) /2. 0; 

PT := FI (T,DF1 ,DF2) ; 

IF PT < P THEN 
ENDL ;= T 
ELSE 

ENDR := T; 

END; 

END; 



. 0) THEN 



. 0 ))* 
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(ENDL ♦ ENDR) /2.0; 



IN VF : = 

END 
ELSE 

IN VF := r: 
END; (* ELSE *) 
END; (* I NVF *) 
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APPENDIX B 



PROCEDURE CON VERTS TRIN3 ; 

( * # ******************** * * ************ lie 



***************** * ******** * ** ********* * 



LABEL 1; 

VAR 

PLACEINT, SIGN ,1 : INTEGER; 

BFDEC : BOOLEAN; 

ACC, TACC, PLACE : REAL; 

PRE VCH , CH ; CHAR; 

BEGIN 

DATA ;= CONCAT (DATA,' ')? 

PRE VCH ; 

BFDEC := TROE; 

INT := 0: 

RL ;= 0.0: 

ACC := 0.0; 

PLACE := 1; 

SIGN := 1; 

FOR I := 1 TO LENGTH (DATA) DO 
BEGIN 

CH •= DATA I * 

IF CH = PREVCH’THEN 
GOTO 1; 

PREVCH := •/' ; 

IF CH IN C' 0* .. '9' ,' . S'] THEN 

CASE CH OF 

*0' , ' 1* 2 ' 3 * , * 4 ' ,'5' , ' 6 ' ,'7*, '8* ,'9*: 



BEGIN 

IF BFDEC THEN 
BEGIN 

RL := RL * PLACE + (OSD (CH)-48) ; 

PLACE := 10; 

END 

ELSE 

BEGIN 

PLACE : = PLACE * 0.1: 

RL .*= RL ♦ PLACE * (0RD(CH)-48) 



RL := RL * SIGN; 

IF A3S (R GER 



* 

* 

* 

* 



PROCEDURE TO CONVERT A STRING 
VARIABLE TO A N333ER 



* 

* 

* 

* 



i • • 

BEGIN 

BFDEC := FALSE 
PLACE := 1 



END; 

END; 

I • * 



END; 

• • 



I • 

BEGIN 



ELSE 

INT 



INT 




END 



1 



# # 



BEGIN 

SIGN := -1: 
PRE7CH := 1 ' 
END; 

' , 1 : GOTO 1; 

END; (*CASE *) 

1 • END * 

END; (* CONVE RISTRI NG *) 
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